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ABSTRACT 


Cerenkov  radiation  is  calculated  for  electron  beams 
which  exceed  the  velocity  of  radiation  in  a  nondispersive 
dielectric  medium.  The  electron  beam  is  assumed  to  be 
bunched  as  emitted  from  a  travelling  wave  accelerator,  and 
the  emission  region  is  assumed  to  be  finite.  Predictions 
include  (a)  emission  of  harmonics  of  the  bunch  rate,  (b)  co¬ 
herence  of  radiation  at  low  frequencies,  (c)  smearing  of 
the  emission  angle  for  finite  emission  regions,  (d)  explicit 
evaluation  of  power  spectrum  in  terms  of  bunch  dimensions. 
The  results  of  theory  are  applied  to  microwave  emission  from 
fast  electrons  in  air  for  different  lengths  of  the  air  path. 
Problems  encountered  during  the  experiment  are  discussed 
and  suggestions  for  continuing  work  are  provided. 
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I .  INTRODUCTION 


Cerenkov  radiation  is  generated  by  a  charged  particle 
moving  at  greater  than  light  speed  in  a  particular  medium. 
Because  the  distribution  of  intensity  of  Cerenkov  radiation 
is  proportional  to  the  frequency,  the  radiated  power  at 
microwave  frequencies  would  be  low  unless  beams  are  intense 
and  bunched  so  that  coherent  radiation  by  many  electrons 
contributes. 

Recent  developments  of  electron  accelerators  for  appli¬ 
cations  such  as  free  electron  lasers  (FED  have  aimed  toward 
high  peak  currents  in  bunches  in  contrast  to  nuclear  and 
particle  physics  applications,  where  low  peak  but  high 
average  currents  are  desirable  to  avoid  saturating  detectors. 
The  high  peak  currents  in  the  new  accelerators  should  yield 
enhanced  Cerenkov  radiation. 

The  intent  of  this  experiment  is  to  find  Cerenkov  radia¬ 
tion  intensity  and  angle  of  emission  in  X  Band  and  K  Band 
regions . 


II.  THEORY 


A.  CALCULATION  OF  THE  POYNTING  VECTOR 

In  the  following  derivaMon,  we  consider  the  Cerenkov 
radiation  produced  in  a  dispersionless  medium  such  as  a  gas 
or  other  dielectric,  by  a  series  of  pulses  of  electrons  such 
as  are  produced  by  a  traveling  wave  electron  accelerator 
(LINAC) .  The  pulses  or  bunches  are  periodic,  che  total  emis¬ 
sion  region  is  finite  and  the  bunches  have  a  finite  size. 

In  determining  the  radiated  power,  procedure  is  to  calcu¬ 
late  the  Poynting  vector  from  fields  which  are  in  turn  ob¬ 
tained  from  solutions  of  the  wave  equations  for  the  potentials. 
Since  the  current  and  charge  densities  entering  into  the  wave 
equations  are  expressed  in  Fourier  form;  the  resulting  fields 
and  radiated  power  also  have  Fourier  components.  In  the  deri¬ 
vation  r  is  the  coordinate  at  which  the  fields  will  be  calcu¬ 
lated,  r'  is  the  coordinate  of  an  element  of  the  charge  which 
produces  the  fields  and  n  is  a  unit  vector  in  the  direction 
of  r. 

We  assume  that  E(r,t)  and  B(r,t)  have  been  expanded  in  a 
Fourier  series,  appropriate  for  the  case  where  the  source 
current  is  periodic.  Then  we  have 

E(r,t)  =  I  E(r,u))  (1) 

and  a  corresponding  expansion  for  5,  where  uj  is  a  discrete 
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variable  and  E  and  B  are  Fourier  series  coefficients.  The 
poynting  vector  S  is  given  by 

^  =  ilx'S  (2) 

u 

and  it  is  easy  to  show  that  the  average  of  S  in  a  direction 
given  by  a  normal  vector  n  is 

i  [  n.?  dt  =  i  y  n.E(r,u;)  X  (3) 

0  ^  '^=-oo 

where  T  is  an  integer  multiple  of  the  period  of  the  periodic 
current . 

-1  '2 

Letting  c  =  (ue)  '  be  the  velocity  of  light  in  the 
medium,  the  wave  equations  for  A,  'j  and  their  solutions  are, 

(7^  -  4-  ^)A(r,t)  =  uJ(r,t) 

3t 

(4) 

2 

(7^  -  ^  4>(r,t)  =  -p(?,t) 

c  9t 

A(r,t)  *  u////D(r-r' ,t-t')3(r' ,t')d^r'dt' 

(5) 

(|)(r,t)  =  i////D(r-r  '  ,t-t ’)  p(r’ ,t ')d^r'dt' 

where  the  Green's  function  D  is  given  by 
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D(r,t)  =■■  ^  6(t  -  r/c 


(6) 


The  vector  potential  A(r,t)  also  can  be  developed  in  a  Fourier 
series  expansion  of  a  form  similar  to  (1)  with  an  expression 
for  the  Fourier  series  coefficients  given  by  (Appendix  B) : 


A(r,cjl 


dt  A  ( r ,  t )  e 


iuj  t 


(7) 


d-r'  5(r',.j) 


2r ”2r  * 


r-r  ’ 


Now  if  we  assume  that  the  observer  is  far  from  the  source 
30  that  r  r',  for  regions  where  the  integrand  in  (7) 

is  important  we  can  let  r-r*'  =  r  -  n*r'  in  the  exponential 
and  r-r'  =  r  in  the  r-r'  factor  in  (7)  obtaining  iwhere 
n  =  r/r) 


A(r ,u) 


u  iwr/c 
4TTr  ® 


J  ) 


- i— n • r ' 

J(r',.^)e 


(8) 


The  Fourier  series  coefficients  of  the  fields  are  obtained 
from  those  for  the  vector  potential  (8)  through  the  usual 
relations  §  =  7  xA  and  E  =  Under  the  conditions 

O  U 

leading  to  (8)  the  field  Fourier  coefficients  are 


§(r,uj)  =  ^nxA(r,u)  (9) 

E(r,  ’>  -  -on  ''R(r.  >  flO) 
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The  poynting  vector  can  now  be  found  by  using  (9)  and 
(10)  in  expansions  like  (1)  cuid  then  substituting  in  (2)  . 
However  it  is  more  convenient  to  deal  with  the  frequency 
components  of  the  radiated  power  by  substituting  (9)  and 
(10)  into  the  expression  for  the  average  radiated  power  (3) . 

ST  /  n-Sdt  =  -  y  ^|nxA(r,j)i  (11) 


B.  FOURIER  COMPONENTS  OF  THE  CURRENT 

The  expression  (7)  for  the  Fourier  components  of  the 
vector  potential  contains  the  Fourier  components  of  the  curren 
density.  Consequently  it  is  necessary  to  examine  the  form 
of  the  current  and  its  Fourier  development.  Assume  the  curren 
is  in  the  z  direction  and  periodic.  If  the  electrons  move 
with  velocity  v,  and  we  ignore  for  the  moment  the  x  and  y 
variables,  the  charge  or  current  functions  should  have  the 
general  form 


f  (z,t) 


=  I  e  ^  I  f(k 


(12) 


Under  the  condition  of  rigid  motion, 


f(z,t)  =  fQ(z-vt) 
It  is  easy  to  show  that 


(13) 


13 


(14) 


f 


2 


where 


0 

Thus  the  restrictions  of  equation  (13)  reduce  the  two 
dimensional  Fourier  series  of  ;i2)  zo  essentially  a  one 
dimensional  series  (14). 

With  (14)  in  mind,  the  current  density  associated  with 
the  electron  beam  from  a  linear  accelerator  should  be  periodic 
in  both  2,  t,  with  a  Fourier  series  expansion  but  the  x  and 
y  dependence  should  be  represented  by  a  Fourier  integral 
form; 


J-(r,t)  =  vo(r,t) 

2 


V 

2z 


:  dk  ~  I  dk 

X  2:'  V 

— OC  —00  * 


r  ^i (k*  r-w • t)  , 

V  "  "0 

K  =—00 

2 


(16) 


where  the  Fourier  components  of  the  charge  density  are 


PQ(r)  =  /  dy  /  /  dze 


-ik -r , 


(r) 


(17) 


—00  —00 


Pg(r)  is  p(r,t)  evaluated  at  t  =  0  and  J  is  assumed  to  be  in 
the  2  direction.  Note  in  (16)  that  k^  and  jj  are  both  dis¬ 
crete  and  from  (14)  at  =  k^v. 
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C.  VECTOR  POTENTIAL 

The  results  of  the  previous  section  can  be  applied  to 
the  evaluation  of  the  vector  potential  and  in  turn  to  the 
fields . 

Let  the  infinite  periodic  pulse  train  be  made  finite,  ex¬ 
tending  from  z  =  -Z '  to  z  =  +Z ' ,  so  that  2Z '  is  the  length 
of  the  path  which  generates  radiation.  Let  6  be  the  angle 
between  n  and  A.  Then  the  cross  product  in  (11;  can  be 
written 


I  n  X  S(r,u)) 


j  ^lijr/c 

=  sin  -  e 

4-n-r 


But 


f  dx '  /  dy '  /  dz  '  e 

-2  ' 


-i-n-r' 


—00  —00 


27  X  •'  V 

•-00  —00  ^ 


=  —00  Z 


0) 


(18) 


f  dx'  /  dy '  /  dz  ' e 


Z'  i?'-(lc-n^) 


—CO  —00 


-z 


=  (2TT)^5(k^-n^  |)5(ky-ny  |)I(Z-) 


(19) 


where 


I(Z  ') 


Z  ' 

/  dz ' 


^  sin  Gz  ' 


(20) 


and 


G 


.  OJ 

K  -  n  — 

Z  2  C 


—  -  n  — 
V  z  c 


And  thus  the  cross  product  term  is 


n  .-<  A(r  ,uj) 


=  sin  6 


(21) 


Note  that  oj  is  a  discrete  variable  but  from  (19)  ,  the 

continuous  variables  k  and  k  become  evaluated  at  discrete 

X  y 

points . 

Returning  to  (17) ,  a  more  symmetric  form  may  be  obtained 
by  assuming  that  CQ(r),  which  is  periodic  in  z  with  period 
Z,  is  actually  zero  between  the  pulses.  Denoting  by  vQ(r) 
the  charge  density  of  a  single  pulse,  which  is  zero  for 
z  <  0  and  z  >  Z  the  integral  on  z  can  be  written 

Z  -ik  z  Z  -ik  z 

/  dze  ^  =  /  ‘^^le  ^  oXd) 

0  ^0  ^ 

<»  -ik  2 

=  f  dze  ^  p^(r)  (22) 

—00 

Then  (17) ,  the  Fourier  coefficient  of  the  charge  density, 
becomes 
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(23) 


Po(k) 


k  ill 


a  re 


pA(r)  = 


where  PQ{k)  is  the  three  diraensional  Fourier  transform  of  the 
single  pulse  described  by  pQ(r).  Substituting  those  expres¬ 
sions  into  (21)  gives  a  final  simple  result  for  the  cross 
product  form: 


n  X  Air  ,ui) 


sin 


^  i^r/c 
4t:  ® 


(k) I  (Z  ') 


(24) 


where 


1(2  •) 


2 

G 


G2 


k  =  (n^-w/c,n^-w/c,^) 

The  components  of  the  Cerenkov  E  and  B  fields  may  now  be 
found  by  substituting  (24)  in  (9)  and  (10) . 

D.  RADIATED  POWER 

The  frequency  components  of  the  average  radiated  power 

are  obtained  by  substituting  (24)  into  (11) .  The  negative 

frequency  terms  equal  the  corresponding  positive  frequency 

terms,  yielding  a  factor  of  2  when  the  summation  range  is 

2 

changed.  Multiplying  by  r  converts  to  average  power  per 

unii-  io'  i.'.;  .;,!:glc,  dl’/'.:..,  v 
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f  0^  n-sdt  =  ^  ^  ^|n  X  A(r,u))  ! 

00 

^  /V 

=  I  W{aj,n)  ,  (26) 

0 

A 

where  W(a),n)  is  defined  to  be 

2  2 

W(uj,n)  =  — ^  sin^e  ^1  o' (k)  j  I^(Z')  (27) 

(4o) ^ 

W(u.',n)  is  the  power  per  unit  solid  angle  radiated  at  the 
frequency  'jj  ,  which  is  a  harmonic  of  the  basic  angular  fre¬ 
quency  Jq  of  the  periodic  pulse  train. 

To  find  P  ,  the  total  cower  radiated  at  the  frequency 
,  W  is  multiplied  by  dfi  and  integrated  over  solid  angle. 
Note  that  =  cos  r,  and  as  9  varies,  G  changes  according 

to  (25)  ,  with  dG  =  -  ^  dn  so  that 

o  z 

dfi  =  d({)  -  dG  (28) 

u 

Nothing  that  the  integral  over  (})  yields  2Tr ,  we  find  the 
result  for  the  total  radiated  power  at  the  frequency  oj  for 
all  angles. 

2  2  G  *’ 

P  =  /  sin^Sio*  (k)  1  ^I^(z)^G  (29) 

u)  4it  c  2^  G* 

The  remaining  integral  over  G  may  now  be  examined.  This 

2 

sin  9  and  factors  may  often  be  slowly  varying  compared  to 


dP 

d^ 
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2 

the  I  (Z')  factor,  which  is  shown  in  Figure  1  for  f  =  8.57 
GHz  land  in  Figure  2]  for  f  =  14.28  GHz].  These  two  figures 
are  the  results  of  (25)  which  was  solved  by  a  Fortran  computer 
program  [Appendix  (D.l)].  For  large  Z',  the  peak  in  I  (Z’) 
becomes  narrow,  and  if  the  integrand  may  be  neglected  outside 
the  physical  range  G'  <  G  <  G" 


I'C 


')  dG 


G' 


4Z 


.2 


,sin“GZ 
'  GZ  ‘ 


■)  dG 


4n-2’  (30) 


Then,  evaluating  the  sin  3  factor  and  PQ{k)  at  the  point 
corresponding  to  G  =  0  (which  is  cos  6  =  =  ^)  shows  that 

3  at  the  peak  of  1(2')  is  the  usual  Cerenkov  angle  3^.  We 
thus  obtain  for  large  Z ' 


? 


7 


sin  3 


c 


2 


31) 


Now  let  22 '/Z  =  ratio  of  the  interaction  length  to  pulse 


spacing  =  N,  the  number  of  pulses.  Also  Z  = 


2t;v 


so  that, 


in  the  large  Z'  limit. 


=  JL 


Wq  V  sin  e^l  p^(k) 


N 


(32) 


To  compare  with  usual  formulations  (32)  is 
to  obtain  the  energy  loss  per  unit  path  length 


divided  by  Nv 
per  pulse: 


dE 


'  ■'3) 
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If  the  pulse  is  in  fact  a  point  charge,  the  Fourier  trans¬ 
form  Oq  (ic)  reduces  to  q  ,  the  total  charge  per  pulse  and 
(33)  is  very  similar  to  the  usual  Cerenkov  energy  loss  formula, 
where  for  a  single  charge  q,  the  radiation  is  continuous  and 
the  factor  ojwq  in  (33)  is  replaced  by  udw .  In  the  present 
case  the  pulse  train  is  periodic  at  angular  frequency  uJq 
and  the  radiation  is  emitted  at  the  harmonic  frequencies 
denoted  by  iy . 

E.  RESULTS  OF  THEORY 

Equation  (29)  and  the  approximate  evaluation  expressed 
as  (32)  form  the  main  result.  Some  consequences  will  now 
be  noted,  and  some  calculations  can  be  made. 

I .  Effect  of  Pulse  Size 

The  spatial  distribution  of  the  charge  in  the  pulse 

appears  in  .'^(k),  which  is  the  Fourier  transform  of  the 

charge  distribution.  One  can  easily  see  from  (25)  and  Figure 
2 

1,  the  peak  of  I  (  Z' )  occurs  at  G  =  0  or  n^  =  c/v.  Thus  at 
the  peak,  iaj/v  =  n  ^  so  that  k,  the  argument  of  pl(k)  ,  is 

Z  w  u 

evaluated  at 


k  =  n  I  (M) 

-¥■ 

We  may  also  define  a  charge  form  factor  F(k,)  as 

p^(k)  =  qF(k)  (35) 
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The  form  factor  F(k)  is  identically  one  for  a  point 


charge,  and  for  a  finite  distribution  F(k)  =1  for  k  =  0. 

Furthermore  F(k)  must  fall  off  as  a  function  of  k 
near  the  origin  if  all  the  charge  has  the  same  sign.  If  the 
pulse  were  spherically  symmetric,  F(k)  would  behave  as  elas¬ 
tic  electron  scattering  form  factors  defined  for  nuclear 

charge  distributions.  In  that  case,  the  mean  square 
2 

radius  <r  >  of  the  charge  distribution  is  given  by  tne  be- 
havior  of  F(k)  near  the  origin. 

2  k^ 

2ZF(k)  -►  i  -  <r  >  -g-  (spherical  pulse)  >'36) 

2 .  Smearing  of  the  Cerenkov  Angle 

For  a  finite  region  over  which  emission  is  allowed, 

2 

namely  if  22'  is  finite,  the  function  I  (z'),  appearing  in 

the  integral  in  (29) ,  will  have  a  finite  width.  Since  the 

2 

peak  height  is  4z '  and  the  area  is  (30) ,  we  can  assign  an 
effective  width  2r  =  area/height  =  tt/z  '  or 

r  =  —  (37) 

2Z  ’ 

Thus  the  radiation  is  emitted  mainly  near  G  =  0 

(which  corresponds  to  9  =  9^)  but  in  a  range  AG  =  tV  .  But 

from  (25) ,  AG  *  ^  An  =  —  A(cos  9)  so  that  there  is  a  range 
c  z  c 


in  cos  9  over  v^ich  emission  occurs: 


Note  that  the  finite  angular  width  of  the  Cerenkov 
cone  angle  in  (38)  has  the  factor  l/uj,  indicating  that  the 
higher  harmonics  are  emitted  in  a  sharper  cone. 

This  result  is  shown  in  Table  1.  Values  are  calcu¬ 
lated  by  use  of  the  seane  Fortran  computer  program  as  in 
Appendix  (Dl) . 


TABLE  1 


Frequency  (GHz) 

O 

O 

I'^min 

1 (cos  5 

8.57 

1.509 

16.129 

0.039 

14.28 

1.509 

12.259 

0.022 

19.99 

1.509 

10.539 

0.016 

25.7 

1.509 

9.249 

0.013 

2  3.36 

1.509 

8.819 

0.011 

37.128 

1.509 

7.959 

0.009 

42.84 

1.509 

7.09 

0.007 

C.  BEHAVIOR  AT  HIGH  FREQUENCIES  RELATED  TO  PULSE  P.ARAMETERS 
To  be  specific  let  the  charge  distribution  for  a  single 
pulse  be  given  by  Gaussian  functions 

2  2  2 

pMr)  =  A  exp( — j  ~  ^ 

a  a  b 


Then  F(k)  may  be  found  as 
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F(lc) 


exp(- 


u2v2 

k  b 
z 


(40) 


Beam  pulse  parameters  could  then  be  determine  by  measuring 
the  Cerenkov  radiation.  For  example,  fast  electrons  from 
an  accelerator  in  air  will  emit  with  a  9^  of  several  degrees 
in  which  case  k^  and  k^  in  (40)  can  be  neglected,  giving 

2  2 

v/b  2  ,2 

F(K)  :  exp( - :  exp(-  ^  (41) 

V 

The  average  power  emitted  into  a  unit  solid  angle  (dP/dil) 
as  a  function  of  angle  (9)  is  shown  in  Figure  3  for  f  =  3.57 
GHz  [and  in  Figure  4  for  f  =  14.28  GHz].  These  two  figures 
are  the  results  of  (26)  which  was  solved  by  a  second  Fortran 
computer  program  [Appendix  (D.2)J. 

Total  power  per  unit  solid  angle  as  a  function  of  angle 
(a)  can  be  calculated  by  a  third  Fortran  computer  program 
[Appendix  (D.3)]  for  many  harmonic  components  in  Figures  5 
and  6.  Figure  5  covers  harmonics  from  3  to  10,  and  Figure  6 
covers  harmonics  from  5  to  10.  For  these  two  figures,  the 
m.aximum  power  and  the  corresponding  Cerenkov  angle  are  shown 
in  Table  2. 


TABLE  2 


^1 

^2 

1  (dP/dn) 

9c 

3 

10 

0.65  W 

5“  .99 

5 

10 

.  5  5 

.1  ^  u 
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Total  Jl’/du 


.T3 


-002  OQC  002  :c^ 


0  (10^  degrees 


Figure  5.  Total  dP/dfi  as  a  Function  of  e  (N  =  3,10) 
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Totdl  UO  W/steradiaii) 
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a 


!  1  f 

I  ‘  • 


c 
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a 


degrees) 


Figure  6.  Total  dP/dfi  as  a  Function  of  8  (N  =  5,10) 
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The  behavior  of  dP/dfi  as  a  function  of  frequency  is  shown 
in  Figure  7  (which  is  plotted  with  the  help  of  the  Fortran 
program  in  Appendix  (D.2)).  The  linear  rise  at  low  frequen¬ 
cies  followed  by  a  fall  off  at  higher  frequencies,  the  peak 
occurring  at 


m 


V 

b 


(42) 


After  inserting  values  of  v  and  b  into  (42)  (Appendix  (C) ) , 

9 

is  found  to  be  126.86  x  10  rad/sec  which  corresponds  to 
in 

the  same  value  as  in  Figure  7. 

The  Cerenkov  angle  is  also  plotted  as  a  function  of  ^ 

in  Figure  8  by  use  of  the  same  Fortran  program. 

Furthermore,  a  different  behavior  would  be  expected  at 

very  high  frequencies.  The  formulation  from  the  beginning 

represents  coherent  radiation  from  all  charges,  not  only  from 

one  pulse,  but  coherence  from  pulse  to  pulse.  F(k)  then 

describes  interference  of  radiation  emitted  from  different 

parts  of  the  pulse,  but  note  that  expressions  (29)  and  (32) 

2  2  2 

will  still  be  proportional  to  q  =  n  e  where  n  is  the  number 
of  electrons  in  a  pulse.  Thus  the  n  dependence  of  indi¬ 
cates  coherence.  But  above  some  high  frequency  such  that 
u),/c  =  2-n/i,  where  i  is  the  mean  spacing  of  electrons  in  the 
cloud,  the  radiation  should  switch  over  to  incoherent  radiation 
from  each  charge  and  P^  should  be  proportional  to  n.  The 
incoherent  radiation  should  then  rise  again  as  a  function  of 


30 


180 


310NU  ai“10S^d3M0d 


(3 

71 

3 


TS 

\ 

a. 

r3 


a 

lu 

3 

O’ 


31 


HRRMONIC  NUMBCR 


III.  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 


It  was  decided  to  find  the  frequency  components  of  the 
average  radiated  power  in  the  X  Band  (8.2-12.4  GHz)  region. 

An  experiment  was  designed  to  measure  average  radiated  power 
as  a  function  of  angle. 

The  equipment  arrangement  is  shown  in  Figure  9.  The 
reflector  used  was  a  thin  sheet  of  aluminum  which  would  allow 
the  passage  of  the  electrons  and  reflect  any  EM  radiation. 

X  Band  antenna,  wave  guide,  detector  (HP  X424A)  are  mounted 
together.  They  were  horizontally  movable  on  this  mount. 
Radiation  was  .measured  using  a  crystal  detector  connected  rc 
an  oscilloscope  (Tektronix  475  A,  200  MHz) . 

Experiments  were  made  for  air  oaths  (L, )  =  66  cm  and  89  cm, 
and  repeated  in  different  days  for  many  times.  Results  cha.nged 
day  by  day  depending  upon  beam  current.  When  the  beam,  cur¬ 
rent  increased,  the  observed  signal  also  increased.  Figures 

10  and  11  show  the  experimental  results  for  Beam  current 
_  8 

equals  2  X  10  Amps. 

Cerenkov  angles  were  calculated  by  help  of  Figures  9  and 
10  (or  Figures  9  and  11)  . 

From  Figure  9,  one  can  show  that 


tan  0^ 


_ X _ 

2  (L^  +  Lj_/2) 


(43) 


where ; 


R£rL£croK 


Figure  10.  Signal  Amplitude  ay  a  luncti' 
(X  Band,  L,  =  09  cm) 


(A«u)  39Hi”10A 


m 


X  =  Distance  between  two  maximum  signal  positions; 

L,  =  Distance  between  LINAC  and  reflector  (air 
path  length) ; 

L2  =  Distance  between  reflector  and  antenna. 

Later  experiments  were  conducted  in  the  K  Band  region  for 
the  same  air  path  lengths  (L^  =  66  and  89  cm)  .  Figures  12 

—  8 

and  13  show  the  experiment  results  for  Beam  current  =2x10 
Cerenkov  angles  were  calculated  for  these  four  experiments. 

The  results  are  shown  in  Table  3. 


TABLE  3 


e'  (XB) 

e;  (KB) 

Lj^  (cm) 

T 

“2 

8,21 

5.94 

89 

9  9 

9.13 

6.97 

56 

So 

To  check  the  correctness  of  the  measurements,  another 
experiment  was  set  up.  The  equipment  arrangement  was  the 
same  as  before  except  this  time  the  reflector  was  rotatable 
(a  different,  bigger  aluminum  thin  sheet  was  used)  and  the 
antenna  was  at  a  fixed  position.  The  experiment  was  repeated 
for  two  different  air  path  lengths.  Rotation  angles  (a)  of  three 
reflector  versus  signal  voltages  are  plotted  in  Figures  14, 

15,  16  and  17. 
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RNGLE  RLPHR 


RNGLE  BLPHR  (Degrees) 


The  relationship  between  9  and  9^  can  be  found  from 
Figure  18.  For  small  values  of  2a  and  26^ 
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c 


2a 


^  rad 
^1 


(44) 


(45) 


From  (44)  and  (4  5)  , 


=  a 


^2  ^^1/2^ 


(46) 


where  a  is  the  rotation  angle  of  reflector  between  two 
maximum  signal  positions. 

Cerenkov  angles  were  calculated  for  various  values. 
Results  are  shown  in  Table  4 . 


TABLE  4 


(cm) 

(cm) 

O 

a 

i°  IKBl 

9°  (KB) 

66 

86 

N  O 

Signal 

Observed 

89 

99 

12 

8.28 

6.04 

109 

114 

11 

7.44 

5.35 
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IV.  DISCUSSION 


A.  EFFECTS  OF  REFLECTOR  AND  ANTENNA 

In  all  experiments,  the  same  X  Band  detector  was  used. 
Table  5  shows  experimental  values  of  6^  and  dP/dfi  assuming 
the  X  Band  detector  is  performing  in  its  Square  Law  region 
where  the  sensitivity  is  400  u'VuW  [Ref.  1]. 


TABLE  5 

2c _ voltage  (mV) _ power  (mW) 


Band 

=  66 

cm) 

9.13 

1.5 

3.75 

X  10“^ 

Band 

(L, 

=  89 

cm) 

8.208 

2 

5 

X 

O 

I 

Band 

=  66 

cm) 

6.67 

0.8 

-  3 

Band 

(Li 

=  89 

cm) 

5.94 

1.4 

3 . 5 

X  10  ^ 

To  make  a  comparison  between  experimental  and  theoretical 
results,  L^  =  89  cm  and  8.2  were  chosen  from  experimental 
data;  and  for  these  values  I  dP/dii  was  calculated  using  the 
second  Fortran  program  (Appendix  C.2).  The  results  are  shown 
in  Table  6 . 

Although  antenna  gains  weren't  included  in  calculations 
of  theoretical  values,  as  seen  from  Tables  5  and  6,  experi¬ 
mental  values  are  much  smaller  than  theoretical  ones.  The 
reasons  for  this  can  be: 

a.  Aluminum  thin  sheet  is  not  an  excellent  reflector. 


TABLE  6 


^1 

dP/dfi  (mW/steradian) 

N=sN^ 

I 

N=N^ 

P  (mW) 

3 

4 

.135 

9.45 

lO'^ 

3 

5 

.223 

15.61 

X  10“^ 

3 

6 

.  30  8 

21.56 

X  lO"^ 

25  .  39 

,  ~-2 

3 

7 

'  VJ 

3 

3 

.424 

29.68 

<  10“^ 

3 

9 

.451 

31.57 

X  10*^ 

3 

10 

.461 

32.27 

X  10“^ 

b.  Antenna  responsivity  is  decreasing  when  distance  between 
radiation  source  and  antenna  increased  (L  =  L]_  +  L2'  - 

To  check  these  conditions,  another  experiment  was  con¬ 
ducted  for  various  L  values  with  a  setup  as  shown  in  Figure 
19  . 

The  results  are  shown  in  Table  7. 


TABLE  7 

L  (cm) 

voltage  (mV) 

power  (mW) 

X 

Band 

80 

130 

3.25  X io~^ 

X 

Band 

160 

20 

0.5  X 10*^ 

K 

Band 

74 

20 

0.5  X  10"^ 

K 

Band 

151 

18 

0.45  X  10“^ 
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In  this  experiment  some  part  of  the  observed  signal  may 

be  caused  by  near  field  radiation,  though,  the  results  now 

seem  more  comparable  with  theoretical  results  and  very  close 
10 

to  I  Pj, . 

N*3  ^ 

To  be  cible  to  make  exact  comparisons  between  experimental 
and  theoretical  values,  antenna  gains  were  found  for  these 
harmonics  of  the  LINAC  [Ref.  2]. 

Calculated  values  are  shown  in  Table  8. 

TABLE  8 

Freq.  (GHz) _ Ant.  Gain  (X  Band) _ Ant.  Gain  (K  Band) 


8.57 

11.57 

— 

11.42 

20.31 

— 

14.28 

30.64 

71.28 

17.14 

42.46 

90.34 

19.99 

55.08 

105.55 

22.85 

68.67 

116.66 

25.70 

81.72 

122.22 

28.56 

94.58 

122.98 

For  =  89  cm,  the  result  of  the  experiment  was  '7^  =  8.2  (from 
Table  1) .  To  make  a  comparison  between  experimental  and 
theoretical  values,  theoretical  value  of  average  power 
emitted  into  a  unit  sold  angle  is  calculated  for  some 
harmonics  of  the  LINAC  by  use  of  the  Fortran  program  in 
[Appendxx  =  o'.  2) 
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The  results  are  shown  in  Table  9. 


TABLE  9 


N 

dP/df2 

(mW/steradian) 

Gain„_  dP/dS3 
(mW/^eradian) 

Gain  g  dP/dQ 
(mW/s€eradian) 

3 

5. 76  X  lO'^ 

.67 

— 

4 

5.62  X  lo"^ 

1.14 

— 

5 

3.33  ■  13"^ 

5  .  33 

6 

8.49  X  10"2 

3.60 

7.67 

7 

6.9  X  lo"^ 

3.80 

7.27 

8 

4.7  X  10"^ 

3.23 

5.48 

9 

2.6  X  lo'^ 

2.12 

3.18 

10 

1.03  X 10"^ 

0.97 

1.27 

As 

seen  from  Tables 

7  and  9,  although  the 

experimental 

values 

include  more  than 

one  harmonic,  the  theoretical 

values 

are  bigger  than  the  experimental  ones. 

B.  CONCLUSION 

Although  in  these  experiments  only  a  very  small  signal 
was  observed,  all  results  indicated  that  the  higher  harmonics 
are  emitter  in  a  sharper  cone.  Nevertheless,  some  refinements 
remain  to  be  done.  Some  specific  suggestions  are  given 
below. 

1 .  Antenna  Sensitivity 

Antenna  sensitivity  has  to  be  measured  as  a  function 

■  requency  arid  di-  -.nc’' . 
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2 .  Detector  Sensitivity 

In  the  K  Band  region,  a  K  Band  detector  should  be 
used.  The  sensitivity  of  the  detectors  should  be  tested. 

3.  Coupling  a  Band  Pass  Filter  to  Detector 

In  the  experiments,  the  X  Bcind  and  K  Band  wave  guides 
have  certain  lower  end  cutoff  frequencies,  but  no  upper  end 
cutoff  frequencies.  If  the  experiments  are  conducted  using 


a  band  cass 


rs  a  I 


the  effectiveness  of  measurements  will  be  increased. 

4  .  Noise  Effect 

The  experimental  area  was  very  noisy.  It  affected 
the  observed  signals  too  much.  Noise  can  be  partially  de¬ 
creased  by  putting  a  detector  into  a  thin  aluminum  box,  which 
would  decrease  the  electric  field  noise  cickuc.  However, 


the  magnetic  field  noise  associated  with  the  LINAC  Klystrons 
is  more  difficult  to  shield. 


APPENDIX  A 


DERIVATION  OF  CERENKOV  RADIATION  FOR  A 
SINGLE  PULSE  OF  CHARGE 

Let  the  pulse  be  described  by 

p  '  (r,t)  =  -  vt)  (Al) 


Both  k  and  ^  are  continuous  variables  in  this  case;  v  is 
z 

again  along  the  z  axis.  If  we  expand  in  terms  of  a  four 
dimensional  Fourier  integral, 


(r,t)  =  (1/{2tt)  )e' 


;t-k*  r) 


( K ,  uj )  d  k  du.' 


\  Al ) 


It  mav  be  shown  that  the  condition  uki)  gives: 


'  (k,oj)  =  2:r6  (u;  -  k^v)  (k) 


[A3) 


where  (k)  is  the  three  dimensional  spatial  transform  of 
p'(r,t)  evaluated  at  t  =  0 .  All  the  fields  have  Fourier 
integrals  rather  than  Fourier  series  expansions  a-'d  the 
energy  radiated  per  unit  solid  angle  becomes 


/  dtn-l  =  I  /  du) 


III /d^r 'dt’e 


i|.(ct'-n-r') 


n^J(r',t')  : 


I  0)  (cj  ,n)  d(i) 


(A4) 
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The  integrand  is  a  synunetric  function  of  oj  so  that 


W(w,n)  =  — ^  ////d^r’dfe^  ,t’) 


167r 


where 


167T 


1  u  2  ^  2  2 

-3-.0  (n^v) 


(A5) 


=  ) /j  I'd^r 0  ’  (r*  ,f  ) 


U^cl 


Now  we  may  write  p’(r',t')  in  a  Fourier  integral  represen¬ 
tation  , 

■wr',t')  =  — ^  ;  d^k'd;^  ’  -  ’  (k' ' 

(2-)  ■’ 

(AT) 

Inserting  (A3)  into  (A7)  and  the  result  into  (A6) ,  the  inte¬ 
gral  over  d^k '  involves  only  exponentials  and  yields 
3  3  -» 

(2ir)  0  (k'  -  ojn  /c)  ,  so  that  (A6)  becomes 


=  /dt'/// /d\'da)'e^^“"''^'^^’o^(k’-wn/c) 


(u,'-k;v)p^ 


k') 


Now  the  integral  over  w '  may  be  done;  because  of  the  i 

function,  u) '  is  evaluated  at  k'v. 

z 

,  i(u)-k'v)t'  ^  ^ 

u  =  /dt’///d  k'e  ^  5^(k'-u)n/c)  P^{k') 
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Now  do  the  integrals  over 
appears  in  the  exponential 


k^,  ky  and  k^ ,  noting  that 
but  k'  and  k'  do  not, 

X  y 


M 


/  dte 


icot ' 


-i^t '  n^v/c 

e  PjJ  ((ijn^/c,ujn^/c,cjn^/c) 


This  may  be  written  as 


M 


dt '  e"^ 


{^n/O 


IA8) 


where 


H  =  1  -  n^v/c 


(A9) 


If  we  let  the  time  interval  be  finite,  frcm  -T  to  ^T,  the 
integral  is  easily  done: 


M  =  -  sin  jHT0Q(nuj/c) 


(AlO) 


sin^-jHT,  ,  ,  ^ 

M  =  4T  - ^I0'(naj/c), 


(All) 


This  result,  equation  (All)  may  be  inserted  in  (A5) 
for  oj.  The  factor  n^^v  is  just  sin  1  where  9  is  the  angle 
between  the  radiation  and  the  beam  axis. 


W(co,n) 


— ^  ii-  a)^sin^94T^ 
16u  ^ 


2 

sin  joHT  I 

- 3~' 

(wHT) 


PjJ  (nw/c) 


2 


(A12) 
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r 


W  is  the  energy  radiated  per  unit  solid  angle  per  unit 
angular  frequency,  oj .  To  proceed  to  the  total  energy, 
multiply  by  dTi  (solid  angle)  and  integrate.  But  n^  =  cos  e 
so  that  dfi  may  be  related  to  dK: 

dfi  =  d(cos  0)dr;  =  -  -  dH  du  (A13) 

V 


The  function  in  equation  (A12)  does  not  contain  ;  so  tha 
integration  over  -3  yields  2~.  Thus: 


W(.^  ,n)  dil 


.  2. 
sin  r 


sin^.jHT 
(wKT)  “ 


dH 


(A14) 


2 

The  sin  jjHT/(cHT)“  factor  in  the  integral  is  peaked  at  H  =  0 

which,  by  equation  (A9;  is  at  n^  =  cos  r  =  ^,  or  the  usual 

Cerenkov  angle,  This  function  is  more  strongly  peaked 

about  H  =  0  for  large  values  of  T,  and  in  fact,  for  large  T 

2 

we  may  evaluated  sin  9  and  .7^  at  the  point  corresponding  to 
H  =  0.  Then  the  integral 


/  dx  sin^ ( ax) / (ax) ^  =  t/a 


may  be  used  to  evaluate  equation  (A14) ,  yielding 


//Wdn  =  ^  ^  2T  sin20^1p'(nw/c)  : 


c'  ^0 


(A15) 


The  emission  was  assumed  to  occur  in  a  time  interval  from  -T 

to  *.  /  <_•  O  O ^  1  ‘  1  j'  Ju  -  .  -  I  ■  'V  i  *-  i.  -  .  ,  li  ^  .  X  I  f 
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and  multiplying  by  v  converts  to  emission  per  unit  path 
length.  Thus  we  obtain,  for  the  large  T  limit: 


.  2 

dxajU  =  47*^^  sin  9^1p^(nu)/c)  I  (Al6) 


2 

where  d  E/dojdx  is  the  energy  emitted  per  unit  path  length 
per  unit  angular  frequency  range  jj . 

The  corresponding  expression  for  T  not  large  is 


d'-E 

dxdoj 


d(jj 


t  ,  ,a,T  2 
_  (_) 


H" 


iin  5  i  Pq  (noj/c) 


2 


2 

sin  a;KT 
(wHT)  ^ 


(A17) 

where  H"  and  H'  are  the  value  of  H  corresponding  to  e  =  0 
and  -j  -  ~  respectively. 

Equations  {A16)  and  {All)  then  describe  the  energy 
radiated  per  unit  path  lenth  and  per  unit  angular  frequency 
range.  For  the  non  periodic  (single)  pulse  the  radiation  has 
a  continuous  frequency  spectrxjm.  For  a  point  charge  q, 

Oq (k)  is  identically  q  and  the  usual  Cerenkov  formula  is 
obtained.  Equation  (Al6)  is  quoted  by  Jelly,  but  only  with 
the  form  factor  corresponding  to  a  uniform  line  charge  of 
length  L'  (Ref.  3]. 
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APPENDIX  B 


DERIVATION  OF  EQUATION  7 

Equation  7  is  derived  for  the  case  in  which  J(r,t)  is 
expanded  in  Fourier  series.  Let  the  Fourier  coefficient  for 
A  be  given  by: 

A(r,>;j)  =  —  /  dt  A(r,t)  (Bl) 

^  0 

Assume  that  the  Green's  fimction  solution  for  A(r,t)  is  given 
as : 

A(r,t)  =  u!;  i’d"r'  '  dt  J(r,t)  D(r-r',t-t'j  iB2} 

where 

0(^,0  =  ^5(t-r/c)  (B3) 

Let  the  current  density  be  expanded  in  a  Fourier  series: 

J(r',t')  =  I  J(r’,i^)  (B4) 

u' 

Then  insert  (B2) ,  (B3)  and  (34)  into  (Bl)  to  obtain 
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A(r,u))  = 


=  i  /  dt  8“''  ///  d\'/df  ^  , 

^0  r  -  r' 


(B5) 


5  (t-f -i  r-r*  i/c)  I  J(?',w') 


10 


Do_  /”  dt',  note  that  t'  appears  in  the  o  function  and 


in  e 


-ioj 't  ■ 


The  result  is  t'  is  evaluated  at  t '  =  t-'r-r' 


A{r,uj) 


i  J'  dt  ///d^r’  ^  ^ 


i  ^  I 


Co  the  integral  on  t,  note  that 


(B6) 


i  /  dt  e 


i  {(jj-'ji ' )  t  _ 


•i)  ,  1 


(37) 


Then  do  the  sum  on  uj ' 


A(r,a)) 


^  "  U-  _  T-  «  I  ~ 


i(jj  I  r-r  '  j  /c 


(B8) 


r-r' 


This  proves  the  desired  result,  (B8)  is  equation  7  as 
used  in  the  main  text. 


/c. 
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APPENDIX  C 
CALCULATIONS 

In  an  earlier  experiment  at  the  Naval  Postgraduate 
School  LINAC,  Leslie  J.  Brown  [Ref.  4]  showed  that  for  small 


(Ci) 


where  e  is  the  electron  energy. 


Ae 


Let  linear  accelerator  be  set  for  - —  =  10 


^-2 


so  that 


"0 


;  V  2  ^  10  ^  rad. 


^^0  -2 

^  =  4 .5  X 10  ^  (C3) 

This  means  that  electrons  are  accepted  4.5%  of  each  cycle. 
Our  LINAC  Klystrons  operate  at  2.856  GHz  (X  =  10.5  cm) , 


.045/2.856  GHz  =  16  p.  sec. 


(C4) 


Therefore,  the  LINAC  produces  bunches  that  are  16  pico~ 
seconds  long.  Incidentally,  it  does  so  in  one  microsecond 
burst  at  the  rate  of  60  bursts  per  second. 
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Figure  20  is  a  sketch  of  the  beam  current  for  simplicity  it 


is  assumed  that  the  current  is  constant  over  its  16  psec.  on 
time,  rather  than  having  some  more  complicated  shape. 

In  theory,  it  is  assumed  that  the  charge  distribution 
for  a  single  pulse  be  given  by  Gaussian  functions 

PQ{r)  =  Aexp(-  x^/a^  -  y^/a^  -  z^/b^) 

2  2 

~  Aexp(-^i^  ^_)  (C5} 

V 

so  it  can  be  assumed  that 

2b  =  .4725  cm 

b  =  .23625  cm  (C6) 


The  condition  for  an  electron  bunch  is  shown  in  Figure 

—8 

21  noting  that  =  2  x  10  Amps. 


=  I 


meter 


1 

TUT 


=  3.3  X 10 


-7 


Amp . 


(C7) 


During  1  psec.  pulse  current  (I2)  is 

I-  =  I,  X - - - r  =  5.55  X  lo”^  Amps  (C8) 

^  ^  60  X lO" 

Using  the  relationship  between  I 2  and  q  (total  charge 

60 


s.^ 

i-v  »_j 


Q 

o 

rO  ^ 


I 

4  ! 

.1  ^ 

^  1 

o  < 


:2  -a- 


ure  20,  Beam  Current  au  a  Function  of  Tiine 


beam 


r’ 


q  =  - - - jr  “  10"^^  Coulomb  (C9) 

2.856  X  10^ 

was  found. 

In  the  present  case,  experiments  are  conducted  in  air  for 
electron  energies  I  100  MeV.  Corresponding  values  are  for 
these  experimental  conditions: 


n 

u 

'0 

B 

C 

V 


refractive  index  of  air  =  1.00036 

permeability  constant  of  air  =  1.256  ^  10  ^  Henry/m 

g 

speed  of  radiation  in  vacuum  =  2.998  x  lo  m/sec 

.999987 

*”0  3 

speed  of  radiation  in  air  =  —  =  2.997  x  lo  m/sec 

8 

speed  of  electrons  =  3  Cq  =  2.997S  ■-.ic  m/sec 


The  relationship  between  dfi  and  area  of  an  antenna  is 


df2 


ant 


area _ 


(CIO) 


The  X  Band  antenna  has  dimensions  4.6  cm x  3.15  cm  and  K 
Band  antenna  has  dimensions  7.2  cm x  5.75  cm.  Results  are 
shown  in  Table  10 . 

The  cut-off  frequency  of  the  X  Band  wave  guide  is  6.56 
GHz,  and  cut  off  frequency  of  the  K  Band  wave  guide  is  14.08 
GHz.  So  all  the  calculations  for  the  X  Band  begins  from  the 

':hird  harmonic,  and  for  K  Band,  from  ‘ho  fifth  harmoric. 
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TABLE  10 


X  Band 
X  Band 
X  Band 
K  Band 
K  Band 
K  Band 


,  2> 
Area  (cm  ) 

Lj^  (cm) 

L2  (cm) 

dS2 

14.49 

88.9 

99.06 

7  X  10“'^ 

14.49 

66.04 

10.16 

-4 

8  X  10  ^ 

14.49 

109.22 

114.3 

-4 

5  X  10  ^ 

41.4 

88.9 

99.06 

20  X  10“* 

41.4 

58.04 

101.6 

-4 

23-10 

41.4 

109.22 

114 . 3 

14.5  X 10 
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APPENDIX  D 


FORTRAN  PROGRAM  LISTINGS 

To  solve  the  equations  (25)  and  (26)  ,  three  Fortran 

computer  programs  were  used. 

2 

1.  FORTRAN  1  solves  I  as  a  function  of  9  for  one  particu¬ 
lar  frequency. 

2.  FORTRAN  2  solves  dP/dr4  as  a  function  of  6  for  one 
particular  frecuency. 

y 

3.  FORTRAN  3  solves  T  dP/di2  as  a  function  of  3. 

N=N, 

In  the  first  and  second  programs,  responses  can  be  found 
for  other  harmonics  by  making  a  change  from  N  =  3  to 
N  =  Desined  harmonic  number.  In  the  third  program,  making 
a  change  from  *K  =  3,10  to  K  =  Desired  har.monic  numbers. 

The  corresponding  program  and  theoretical  syrhols 
are  shown  in  Table  11. 


TABLE  11 


Program  Symbol 

^1 

PI 

WO 

V 

C 

Q 

i  ,  1 U  2  4  'J  L)  o  4  4  'j 

-0.01000480877 


Theoretical  Symbol 
Z  ’ 


'^0 

V 

c 

^222 
M  u)q  V  q 

.72  c- 

1  ,‘^ob,  2 
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1)  FORTRAN  program  1 


//AHMETI3  JOB  ( 2075 , 1000 ) t* AT  ILL AX* , CL  ASS* A 
//  5XEC  PRTXCLGP, REGION. G0*300K 
//FORT.SYSTN  00  • 

REALMS  TITLE!  121 
PEAL*8  LABEL/'  •/ 

REAL  X(  llOOlt  YdlOOIfOKllOO) 

REAL  WOtVtCtZl.PI «N«gt 

105  FORMAT! 6Aa) 

Z13.4445 

PI  *3.  14159263 

N*3. 

W0*17. 94477724*! 10**9) 
V*2.997886027*!10**8 ) 

C* 2. 997025892*! 10**8) 

X(  U=0. 

00  10  1=1,991 

OKI  »*X(n-22. 

A*N*^i*?Wo9v-WO*COS!QI/C) 

8l  =  S  IN!  A  i*2. 

G=<N*(WO/V-WO*COS!QJ/C1 ) 

F*B1/G 
Y(I)»F**2 
X!  I*1)  =  X(I )*0.043 
10  CONTINUE 

CALL  DR AW ( 990,01, Y, 0,0, LABEL, title, 
10. ,0.,0.,0.,0.,0.,4,4,1, LA  ST ) 

WRITE<6,25) 

WRITE! 6, 100) ( 01! I) ,Y(I I ,1*1,990,10) 
100  FORMAT!  l2X,ei6.7, 10X,E16.7) 

25  FORMAT!  »1'  ,  15  X  ,  •  AiMGL  =  •  ,  2 IX,  •  I  SOU  AR  E  *  ) 
STOO 
END 

//G0.SYSTN  DO  * 

ANGLE  VERSUS 
I  S3UARS(N*3) 
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2)  FORTRAN  program  2 

//SAGLAM3  JOB  (  2075 . 1000 ) . * AHHETX* ,CLA S5= A 
//  EXEC  FPTXCLGP, REGION. G6=300K 
//FQRT.SYSTN  00  ♦ 

PEAL*8  TITLE! 12) 

REAL*a  LABEL/*  */ 

REAL  X(  1100), YCllOC) ,01(1500) 

REAL  WO,V,C,Z 1,P I,N.Q 
PEAD(5,105) (TITLE( J) ,J=1,6) 
BEA0(5,105)(TITL6( J) ,J*7,12) 

105  F0PMAT(6A8) 

Zl=.4445 
PI*3. 14159263 
N=3. 

W0=179. 4477724*1 10**9) 
V=2.997886027*(10**8 ) 

5YP=2. 713281828 
C=2.997025892*( 10**8) 

X( 1 )=0. 

on  10  1=1,991 

OKI  )  =  X(I)-22. 

Q=Pr/l8C.*Ql ( I ) 

A=N*Z1*<  v*0/V-W0*CCS(Q)/C) 

B1=S  IN( A  ) 

B=61**2 

G  =  (N*Zl*(WO/V-rtO*COS(C)/C)  )**2 
0=2. 1024966449*( ( N*S I N ( Q ) *Z 1 ) **2 ) 
|l|-^.JjJ00480a77*(i N*COS( 0) )**2 ) 

F  =  fl/G 

Y( I )=E*F*0 
X(l+l)«X(I)+0.043 
10  CONTINUE 

CALL  0RAW(990  ,01 , Y ,0 ,0 , LABE L , T I  TLE , 

10. ,0. ,0»,0»,0«,0.,^,4, i,LAST) 

W»  ITS  {ij,25) 

WRITS  (6, 100)  (OK  I  )  ,Y(  I)  ,1=1  ,990,  10) 

IQO  FORMAT) 10X,Elb.7,5X. £lo.7) 

2  5  FORMAT!  *1'  ,I5X,  'A.NGLc*  ,  16X ,  •  OP  /  DCMEGA  •  ) 
STOP 
ENO 

//GO.SYSIN  00  * 

angle  versus 

OP/OOMFGA  (N=3) 
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3)  FORTRAN  program  3 


//S=LMA  JOB  (  2n75,lCC0) f • SELMAX* ,CLASS=A 
//  EXEC  FaTxCLGP.PEGICN.GO*30CK 
//POR'^.SYSIN  00  ♦ 

REAL*8  T!TLE( 121 
RFAL*8  LABEL/*  •/ 

P^AL  X(  1100) 1 Y<1100J ,Q 1 < 1500 ),PZ( 25) 

RPAL  WO»V,C,ntPI,N,Q,T(  1100) 

REAOC  5, 105) (TITLE( J) , J=l .6) 

REA0(  5,105)  (TITLE(  Jl  ,J=7,12) 

105  FORMAT(6Aa) 

Z1=.4445 
00  12  L= 1,055 
T(L) =0. 

12  continue 

PI =3. 14159263 
V«0=179. 4477724 
V=  2.  9978d6027 
EXP*2 .7182318  28 
0=2.997025892 
DO  10  K=3,10 
N=‘:LrAT  (K  ) 

X(  1)=0. 

00  20  1=1,995 
OKI  )=X(  I) 

Q=PI /13G.*Q1( I ) 

A^N^Zl^i WO/V-WO*COSCQ)/C) 

B1=SIN<  A) 

e=81**2 

G=<N*Zl*(hO/V-wO*COS(Q)/C) )**2 
0=2.10249664**9*<  (  N*S  1  N <Q  )  ♦  Z1  )  **2  ) 

E l=-0.0  1000480077*1 (N*CO SI Q) )**2) 

E  =  EXP**E1 
F  =  6/G 

Y  (  I  )*E*F*0 

X( I  +  1)  =  X  (I )+Q.C5 

T( I )=T(  I  )+Y( I ) 

20  CONTINUE 
10  CONTINUE 

CALL  0RAW(995 ,01, T, 0,0, LABEL, TI TLE, 

10. ,0. ,0« ,0.,0.,0.,4,4, 1, LAST ) 

WR1TE(6,25) 

WR  ITE(6, 100)  (OK  I  ) ,  T(  I  )  ,  1  =  1 ,990,  10) 

100  F0RMAT<10X,E16.7,5X,E16.7) 

25  FORMAT! • 1* , 15 X  ,* ANGLE *, 18X, • DPT/ DOMEGA* ) 
STOP 
END 

//GO. SYS  IN  OD  • 

ANGLE  VERSUS 

TOTAL  DP/OOMEGA  <N»3,10J 
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